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We unveil novel spin textures in an itinerant fermion model on a frustrated triangular lattice in the 
limit of low electronic density. Using hybrid Monte Carlo simulations on finite clusters we identify 
two type of nanoscale spin textures in the background of 120° order : (i) a planar ferromagnetic 
cluster, and (ii) and a non-coplanar cluster with spins oriented perpendicular to the 120° plane. 
Both these textures lead to localization of the electronic wavefunctions and are in-turn stabilized 
by the concomitant charge modulations. The non-coplanar spin texture is accompanied by an 
unusual scalar chirality pattern. A well defined electric charge and magnetic moment associated 
with these textures allow for their easy manipulation by external electric and magnetic fields - a 
desirable feature for data storage. We identify a localization-delocalization behavior for electronic 
wavefunctions which is unique to frustrated magnets, and propose a general framework for stabilizing 
similar spin textures in spin-charge coupled systems. 

PACS numbers: 71.27.-l-a, 71.10.-w, 72.15.Rn, 75.10.-b 


Introduction — Extended magnetic objects that are 
large enough to be stable against thermal and quantum 
fluctuations are in common use to densely store infor¬ 
mation. In materials with immobile magnetic bits, such 
as for instance patches of magnetization in hard ferro- 
magnets, there are essentially two ways to address in¬ 
dividual bits. If one wishes steer clear of the overhead 
of wiring each bit individually, the choice that remains 
is to physically move the storage medium with respect 
to the reading/writing device. In materials with mobile 
magnetic information carriers the required mechanically 
moving parts can in principle be avoided by propagat¬ 
ing instead the mesoscopic magnetic bits themselves, for 
instance by allying an electric field in a race-track mem¬ 
ory set-up [U, Q • In tbe past years this concept has given 
a strong impetus to theoretical and experimental research 
on mobile magnetic textures, such as single magnetic 
skyrmions, polarons, bubbles and domain walls |3l-lll|. 
On the theoretical side, the focus so far has been on find¬ 
ing novel spin textures in purely magnetic systems. The 
candidate models are the anisotropic Heisenberg models 
in the presence of external magnetic field 0, H, IH, [ll] ■ 

In this Letter, we show that such novel spin textures 
resembling magnetic polarons can emerge in frustrated 
antiferromagets made out of spins and charges on a tri¬ 
angular lattice, which in their phase diagram are close to 
a region of phase separation. The results are obtained 
via unbiased Monte Carlo simulations on finite lattices, 
and crosschecked with variational calculations. We find 
that ferromagnetic polarons with two different spin ori¬ 
entations are stabilized as a consequence of competing 
interactions on the frustrated lattice geometry. These 
spin textures derive their stability from the fact that they 
appear together with a charge modulation. As opposed 
to most other magnetic skyrmion and domain-wall tex¬ 
tures these spin-textures are thus heavily charged and 
mobile. This implies that they can easily be manipu¬ 


lated with applied electric and magnetic fields, which is 
a useful feature for data storage. From a fundamental 
standpoint these results provide a connection between lo¬ 
calization physics, frustrated magnetism and novel spin 
textures. From the charge profiles we infer the presence 
of an interesting localization to delocalization behavior 
for electronic wavefunctions that is specific to frustrated 
lattices. 

Model and Method - The Hamiltonian that we con¬ 
sider consists of localized classical moments Heisenberg- 
coupled to nearest neighbors and Kondo-coupled to itin¬ 
erant fermions residing on a triangular lattice. The inter¬ 
play of conduction electrons and the frustrated geometry 
give rise to many interesting effects, such as for instance 
a multiple-Q ma gne tic order [T3 - [T^ , the anomalous Hall 
effect [3, Ea, , coupled spin-charge phases , 

partially disordered phases [^, . We work in the limit 

of strong Hondo coupling and focus explicitly on the low 
electronic density limit. In the strong coupling limit the 
electronic spin gets slaved to the local moment and the 
Hondo model reduces to a double-exchange (DE) model. 
The resulting Hamiltonian for spinless fermions is given 

by, 

i? = - X X 

(U> (U> 

where Cj (c|) is the electron annihilation (creation) oper¬ 
ator with spin parallel to the local magnetic moment Si. 
The angular brackets in the summations denote the nn 
pairs of sites on a triangular lattice. Jaf is the strength 
of the AF coupling between the localized spins. The pro¬ 
jected hopping is given by, tij = to[cos(0i/2) cos(0j/2) -|- 
sin(0i/2) sin(0j/2)e“‘('^‘“‘^^b]^ where 9i {(jii) is the polar 
(azimuthal) angle for localized spins at site i, and to is 
the bare fermionic hopping amplitude between nn sites. 
Here onwards we set to = 1 £^s the unit of energy, and 
therefore the only free parameter in the model is Jaf- 
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FIG. 1. (Color online) (a) Total energy per site for different 
values of Jaf- The empty circles are the Monte Carlo re¬ 
sults, the solid line is obtained from variational calculations 
for planar FM cluster (see text), and the star symbols are re¬ 
stricted Monte Carlo simulations where some sites were forced 
to retain the 120° order. The dashed line with positive (neg¬ 
ative) slope represents the energy of a long-range ordered FM 
(120° ) state. Two typical Monte Carlo snapshots showing, 
(b) a planar FM cluster and (c) A non-coplanar FM cluster, 
residing in the 120° background. 


The model has been extensively studied in the context 
of colossal magnetoresistance l26l-[^ , and more recently 
for frustrated itinerant magnets [la l29l - K^ . 


Here we focus on the limit of very low electronic filling 
in order to understand the effect of itinerant electrons 
on a frustrated magnetic state. Note that the DE model 
alone {Jaf = 0) leads to a ferromagnetic ground state 
even for a single itinerant electron. On the other hand, 
in the absence of electrons the ground state supports a 
120° order for all non-zero values of Jaf- While the 
overall ground state phase diagrams for this model sug¬ 
gest a tendency towards electronic phase separation in 
the low-density regime, the detailed real-space nature of 
the states is not known (3^ . In particular, one can ask 
if certain special spin textures are induced as a conse¬ 
quence of a single electron trying to gain kinetic energy 
in an otherwise frustrated magnetic structure. In order 
to find the answer, we make use of a hybrid Monte Carlo 
method which combines the classical Monte Carlo for 
spins with numerical diagonalization for fermions [3^ . 
The method is numerically exact, and has served as a 
very useful tool for exploring the physics of spin-charge 
coupled systems . The simulations begin at high tem¬ 
perature with a random spin configuration, and the the 
temperature is then decreased in small steps [s^. The 
results presented here correspond to low temperatures, 
ksT - 0.002to. 


Results and Discussions: The total energy per site for 
the case of one electron on a 12x12 lattice is shown in 
Fig. in a). As expected, the MC energy (circles) interpo¬ 
lates between the FM and the 120° energy asymptotes 
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FIG. 2. (Color online) (a)-(d) Spin and charge configura¬ 
tions at different Jaf values as obtained from MC simulations 
where spins were restricted to remain in a plane. The color of 
the spin represents the charge density at that site. The charge 
tends to initially localize and then delocalize upon increasing 
Jaf values. 


which are plotted as dashed lines. A careful analysis of 
the low-temperature real-space spin configurations indi¬ 
cated the presence of two interesting spin-textures in the 
intermediate Jaf regime, 0.01 < Jaf < 0.05. The first 
one is a planar FM cluster (pFC) and the second is a 
non-coplanar FM cluster (nFC) with additional meron- 
like texture (see Fig. [T](b)-(c)). These two spin textures 
are almost degenerate in energy over the entire interme¬ 
diate Jaf range. Note that it is difficult to identify the 
detailed nature of such textures from a calculation of spin 
structure factors, as the majority of the background re¬ 
tains the 120° order and therefore the structure factor 
will be dominated by the 120° state. Taking hints from 
the textures obtained in the Monte Carlo simulations, 
we set-up restricted Monte Carlo simulations in order to 
better understand the nature of these spin patterns. 

In the first scheme we restrict the spins to lie in a sin¬ 
gle plane so that the spins are effectively XY type. This 
can also be thought of as a model for easy plane trian¬ 
gular magnet. The ground state spin and charge con¬ 
figurations obtained from these simulations are shown in 
Fig. [21 While the large regions of the lattice still fol¬ 
low the 120° pattern, a cluster of aligned spins emerges 
spontaneously. The charge density is large over these 
aligned spins, and is essentially zero elsewhere. This pro¬ 
vides a simple picture of electronic self-trapping in the 
120° background: a single electron tries to make a FM 
cluster so that it can move freely over this cluster, how¬ 
ever there is a price to be paid in terms of the Jaf- 
Therefore, one expects that the size of the FM cluster 
should decrease upon increasing Jaf- Indeed, this is re¬ 
flected in the real-space plots shown in Fig. |21 However, 
even though the size of the FM cluster decreases mono- 
tonically, the charge patterns indicate a non-monotonic 
localization-delocalization crossover. To investigate fur¬ 
ther the variations of FM cluster size and the charge 
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FIG. 3. (Color online) The radius (red stars), Vc, and the 
angle relative to the 120° state (blue circles), (jic, of the planar 
FC as a function of the AF coupling strength Jaf, as obtained 
from Monte Carlo simulations where spins are restricted to 
be in a single plane. The filled squares represent the charge 
density, no, at the central site of the planar FC. The charge 
density is scaled up by a factor of 5 in order to clarify the 
variations, (b) The variation of the radius, Tc, of the non- 
coplanar FC and no inferred from a second set of restricted 
Monte Carlo simulation (see text). 


localization, we perform variational calculations by al¬ 
lowing for FM clusters with different radii, and different 
orientations w.r.t. the 120° order. This is achieved by 
orienting all spins inside a ring of radius r in same direc¬ 
tion, while the remaining spins retain the 120° pattern. 
For each value of Jaf the FM cluster is varied in ra¬ 
dius as well as orientation to obtain the minimum energy. 
The results are plotted in Fig. [3Ka). While the radius of 
the FM cluster decreases monotonically with increasing 
Jaf, the orientation shows a non-monotonic behavior. 
Tc —>■ oo, (j)^ = Q correspond to a FM phase, and Tc = 0, 
= 0 is the 120° phase. Thus the variational calcula¬ 
tion captures naturally the two limiting phases. In the 
same plot we show the charge at the central site, noj of 
the FM cluster, as a function of Jaf- There is a non¬ 
monotonic variation in ng which indicated an interesting 
localization-delocalization crossover. While the localiza¬ 
tion is gradual with charge at central site increasing to 
its maximum possible value around Jaf ~ 0.05, the de- 
localization is abrupt near Jaf ~ 0.15. 

As mentioned earlier, another spin texture is almost 
degenerate in energy with the pFC. This is a non- 
coplanar spin arrangement where a cluster of spins ori¬ 
ents normal to the 120° plane. To analyze this structure 
and its stability further we perform another set of re¬ 
stricted Monte Carlo simulations. We force the boundary 
layer of spins to retain the 120° order, and a central spin 
to point perpendicular to the 120° plane. The remaining 
spins are updated using the Monte Carlo simulations [s^ . 
The resulting ground state spin configurations are shown 
in Fig. 0] (a)-(b) for two values of Jaf- Here FM clusters 
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FIG. 4. (Color online) (a)-(b) The non-coplaner spin config¬ 
uration at Jaf = 0.005 and Jaf = 0.03 obtained within re¬ 
stricted MC simulations, where the boundary layers are forced 
to retain 120° order and a central spin in fixed to be normal 
to the 120° plane. The magnitude (sign) of Sz{i) is shown 
as the size (color) of the circles, (c)-(d) The scalar chirality 
pattern for the two values of Jaf- 


point perpendicular to 120° spin configuration. On the 
boundary of the FM cluster, ^-components of the spins 
on different rings oscillate in sign before they vanish ap¬ 
proaching the boundary. The z-components of all spins 
at each site are also shown in panels (a) and (b). The size 
of the filled circles represent the values of z-components 
and colors (red or blue) represents their sign. The corre¬ 
sponding charge distributions (not shown) are again cor¬ 
related with the formation of FM clusters. The radius of 
the non-coplanar FC, together with the density at cen¬ 
tral site uq are shown in Fig. [2Kb). Similar to the case of 
planar FC, the cluster size decreases monotonically upon 
increasing Jaf and the variation in uq is non-monotonic 
Note that the nFC texture is not simply a pFC texture 
rotated by 7r/2, but contains interesting modulations of 
the ^-components of the spins (see Fig. 0]). In order to 
highlight the noncoplanar character of the nFC texture, 
we compute the associated chirality patterns. These are 
shown in Fig.0|(c)-(d), where Xijk = Si-(Sj xS^) for each 
triangle with vertices ij and k is shown in the center of 
the triangle. The size of the circle represents the magni¬ 
tude of the chirality and the color indicates the sign. The 
positive and negative chiralities are intertwined to result 
in a zero net chirality. The textures represent examples 
of localized chirality quadrupoles. 

In order to test the stability of these textures in a re¬ 
alistic scenario, it is important to look beyond the sin¬ 
gle electron case. While a genuine finite-density system 
would require us to go to larger lattices which is not 
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FIG. 5. (Color online) The spin-charge textures for the case 
of two electrons for (a) Jaf = 0.03 and (b) Jaf = 0.1. (c) 
Variation of energy with distance d between the centers of the 
two rc = 1 planar FM cluster for Jaf = 0.03. (d) Same as 
(c) for Jaf = 0.1, where only one spin is deviated from the 
120° order. 


feasible within the exact simulation scheme, we present 
results for two electrons on a 12 x 12 lattice [^. Fig.[5Ka) 
and (b) display the ground state spin-charge configura¬ 
tion obtained from the MC simulations with XY spins 
with two electrons for Jaf = 0.03 and Jaf = 0.10, re¬ 
spectively. The spin-charge textures obtained in the one- 
electron case retain their identity and their separation is 
determined by the value of Jaf- It is important to note 
that two electrons do not lead to a larger FM cluster 
with uniform charge density. For a better understanding 
of the two-electron results in terms of single electrons, we 
study the effective interaction between two pFC textures. 
This is achieved by assuming a fixed radius for the two 
pFCs and plotting the total energy as a function of their 
separation d. The resulting plots are shown in Fig. [Sjc) 
and (d) for the two representative values of Jaf- For 
Jaf = 0.03, we find that the energy is minimum at a 
separation of d = 3. Indeed, the separation found in the 
Monte Carlo simulations is close to 3 lattice spacings. 


For Jaf = 0.10, the energy keeps decreasing monotoni- 
cally with increasing d, indicating an effective repulsive 
interaction between two pFC textures. The Monte Carlo 
results are consistent with this simple picture and the 
two textures reside away from each other (see Fig. [5Kb)). 

A remark is in order concerning the generality of such 
spin-charge textures in Kondo-lattice Hamiltonians. A 
thermodynamic limit argument suggests that for any fi¬ 
nite electronic density inside the phase separation regime, 
the system will show macroscopic phase coexistence. 
However, this assumes that the boundary contributions 
are negligible compared to bulk contributions. This as¬ 
sumption does not hold if the lattice is not infinite or if 
the volume fraction of one of the phases is much smaller. 
Therefore, such spin textures can be generally expected 
if the density lies close to one of the phase separation 
boundaries in mesoscopic systems. 

Conclusion: We show that electronic Hamiltonians in¬ 
volving coupled spin and charge degrees of freedom are 
capable of stabilizing novel spin-charge textures in the 
ground state. We study a prototype model for coupled 
spin and charge degrees of freedom using a numerically 
exact Monte Carlo simulation method. In the low-density 
limit, we find that two novel spin-charge textures are sta¬ 
bilized as a consequence of competing interactions on a 
frustrated lattice. These are the planar FM cluster and 
the non-coplanar FM clusters, both of which exist to¬ 
gether with a localized charge-density profile. We also 
propose a general scheme for stabilizing such spin tex¬ 
tures in electronic Hamiltonians involving coupled spin 
and charge variables. The search for stable and mobile 
spin textures in magnetic systems is a rapidly evolving 
research direction due to the potential applications in 
data storage devices with easy read/write operations. 
In contrast to the conventional mechanism of stabiliz¬ 
ing novel spin textures in magnetic models, which in¬ 
volves anisotropic interaction^ longer-range interactions 
and external magnetic fields [3, |^, [13, [l3| ; the electronic 
route presented in this letter can open up new possibili¬ 
ties. 
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In this supplement we include the following material: (i) the details of the hybrid Monte Carlo 
method and the restricted simulation schemes used in this work, (ii) results on the effect of finite 
lattice sizes and that of increasing number of electrons on the stability of the spin textures obtained 
in our study. 


A. Hybrid Monte Carlo Method 


The model Eq. (1) of the main text is investigated us¬ 
ing a hybrid Monte Carlo (MC) method which combines 
the classical Monte Carlo for spins with the fermion diag- 
onalization. The solution of a new Schroedinger equation 
is required at each Monte Carlo update step in order to 
obtain the electronic contribution to the total energy of 
a given classical configurations of localized spins. This is 
achieved by numerical diagonalization of the Hamiltonian 
matrix at each Monte Carlo step. The standard Metropo¬ 
lis algorithm is then followed for spin updates. The com¬ 
putational time with number of lattice sites scales as 
{N from Monte Carlo and from matrix diago¬ 
nalization). Most of the results presented in this work 
are obtained on clusters with N = 12^ sites. Typically 
~ lO'^ MC steps are used for equilibration at each tem¬ 
perature, and a similar number of steps for computing 
thermal averages over classical spin variables. We begin 
the simulations with a random spin configuration and 
at a high temperature ^ 0.3) and reduce stepwise 
in about 30 temperature steps to very low temperature 
{UbT ~ 0.002). The results at these low temperatures 
are taken as representative of the ground state proper¬ 
ties. 

Since the focus of our work is the formation of 
nanoscale textures in the low density limit, the usual 
thermodynamic quantities such as the structure factors, 
etc. are not very useful. Since most of the region retains 
a 120° order, the structure factor will only reflect the 
presence of that order. Therefore, the identification of 
textures requires a careful real-space analysis. In order to 
achieve this we have made use of restricted Monte Carlo 
simulations, where the type of restrictions we impose on 
spin variables is motivated by the structures obtained 
from fully unbiased simulations. From the full simula¬ 
tions we could identify the tendency for formation of two 
type of textures, in the background of the 120° order: 
(i) the planar ferromagnetic cluster, where the spins of 
the cluster align along a common direction within the 
120° plane, and (ii) an out of plane ferromagnetic cluster 
with modulations in Sz, which has a cluster of spins lift¬ 
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FIG. 1. Color-map plot of total energy as a function of polar 
and azimuthal angle of the 7-site FM cluster. The minimum 
energy corresponds to 0 = 7r/2 and (f = n. 


ing out of 120° plane. We used two different restricted 
Monte Carlo schemes in order to further check the ener¬ 
getics of these spin textures. 


B. Restricted Monte Carlo Simulations 

In order to check the existence and stability 
of the planer ferromagnetic cluster (pFC) in the 
120° background, we perform MC simulation by limit¬ 
ing the spins to remain in XY plane. Indeed, the energy 
obtained in these planer simulations is further lowered 
compared to that obtained in the unrestricted simula¬ 
tions (see Fig. 1 in main text). 

As mentioned earlier, Monte Carlo simulations also in¬ 
dicate the stability of textures where a cluster of spins 
gets aligned and points perpendicular to the 120° plane. 
In order to verify the stability of such a non-coplaner 
ferromagnetic cluster (nFC) we perform another set of 
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FIG. 2. Spin and charge configurations at Jaf ~ 0.03 with 
one electron on clusters with, (a) 6^, (b) 9^ and (c) 12^ sites 
respectively. The arrow directions indicate the orientation 
of planar spins and the color on the arrows represents the 
electronic charge density. 


restricted simulations. In these simulations, one of the 
spins is forced to point along z direction and all the spins 
that are located at a distance of r/a > 4 (a being the 
nearest neighbor distance) from the fixed spin are forced 
to retain the 120° order in XY plane. The choice of the 
cut-off of r/a > 4 is motivated by the fact that the unbi¬ 
ased calculations lead to a localization of charge density 
of the electron over sites that are well within the cut-off. 
Therefore, in the ground state 120° pattern is expected to 
be retained in these charge-free regions. The remaining 
spins are then updated using the standard Metropolis al¬ 
gorithm. These simulations also improve the total energy 
of the system compared to unrestricted simulations. In 
fact, the value of total energy of configurations contain¬ 
ing pFC and nFC textures lie very close to each other. 
As discussed in the main text, the nFC texture also has 
modulations in the z component of the spins, which seem 
to crucial for the stability of the texture from an ener¬ 
getic point of view. Therefore, it is important to empha¬ 
size that the nFC is not simply a pFC texture rotated by 
7r/2. This is explicitly confirmed by the color-map plot 
of total energy shown in Fig. SI. The total energy of 
a 7-site cluster is plotted as a function of polar angle 9 
and azimuthal angle (j). The minimum energy is obtained 
for one specific choice of the {9, (j)) pair, indicating an 
anisotropic nature of the pFC. This is in contrast to the 
standard magnetic polarons living in the paramegnetic 
background. 


C. Effect of Lattice Size and Electron Numbers 

In Figl2l we present results on the effect of lattice size 
used in simulations. The emergence of pFC texture is 
already clear for = 6^, N being the number of sites. 
The cluster consists of 7 sites with their spins almost 
aligned. The FM cluster size remains independent of the 
lattice size N, as shown in the figure. Therefore, the 
formation of such texture is accurately captured even on 
small lattices. 

Another very important issue is to verify that such 
individual texture survive when the number of electrons 

(a) J,p=0.03 (b) Lf=0.07 

0.3 



FIG. 3. (a)-(d) Spin and charge configurations on a lattice 
with N = 12^ for one, two, three and four electrons respec¬ 
tively. The values of Jaf are indicated in the plot. The 
arrow directions indicate the orientation of planar spins and 
the color on the arrows represents the electronic charge den¬ 
sity. Note that a larger FM cluster is not formed and these 
textures remain stable as individual dressed particles. 


increases. In order to test this explicitly, we perform 
simulations for different number of electrons on a 12^ 
lattice. The results are shown in Figl31 Clearly, a pFC 
exists as an isolated texture and number of pFC increases 
upon increasing the number of itinerant electrons. It 
is important to note that a larger uniform-density FM 
region hosting all the electrons does not appear as the 
number of electrons is increased. This suggests that at 
low electronic density, where the mean electron-electron 
distance is large, such textures can be realized. 












